X-ray activity phased with planet motion in HD 189733? 
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ABSTRACT 

We report on the follow-up XMM-Newton observation of the planet-hosting star HD 189733 
we obtained in April 2011. We observe a flare just after the secondary transit of the hot Jupiter. 
This event shares the same phase and many of the characteristics of the flare we observed in 
2009. We suggest that a systematic interaction between planet and stellar magnetic fields when 
the planet passes close to active regions on the star can lead to periodic variability phased with 
planetary motion. By mean of high resolution X-ray spectroscopy with RGS we determine that 
the corona of this star is unusually dense. 
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Introduction 
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The evidence of star-planet interaction (SPI) 
in the X-ray band is a lively matter of debate. 
To first order, close-in giant planets (also known 
as "hot Jupiters") should affect their host stars 
through both tidal and magneto-hydrodynamical 
effects (cf. Cuntz et al. 2000 and Ip et al. 2004). 
Both effects should scale as the —3 power distance 
between the bodies (Saar et al. 2004). Kashyap 
et al. (2008) showed that stars with hot Jupiters 
are statistically brighter in X-rays than stars with- 
out hot Jupiters. On average Kashyap et al. (2008) 
observed an excess of X-ray emission by a factor 
of 4 in the hot Jupiter sample. Interplay between 
the magnetic fields of the hot Jupiter and the star 
may be the source of this difference. This could be 
due to interacting winds and magnetic fields or in- 
directly by enhancing the stellar dynamo. On the 
other hand, Poppenhaeger et al. (2010) (and Pop- 
penhaeger & Schmitt 2011) have found no statis- 



tical evidence of X-ray SPI as claimed by Kashyap 
et al. (2008). 

The system of HD 189733 offers a unique en- 
vironment to study SPI effects in X-rays and dis- 
entangle them from proper coronal activity. It is 
composed of a K1.5V type star (at only 19.3 pc 
from Sun), and a M4 companion at 3200 AU from 
the primary, orbiting on a plane perpendicular to 
the line of sight. It hosts a hot Jupiter class planet 
(HD 189733 b) at a distance of only 0.031 AU with 
a orbital period of ^ 2.22d (Bouchy et al. 2005). 

In 2009 we observed the eclipse of the planetary 
companion to HD 189733 with the goal of study- 
ing star planet interactions (SPI) in the case of a 
hot Jupiter (Pillitteri et al. 2010, hereafter Paper 
I). We observed a softening of the spectrum in 
strict correspondence of the planetary eclipse, and 
a flare which occurred 3 ks after the end of the 
eclipse of the planet. The non-detection of the M 
type companion is a strong constraint on the age 
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of the system at > 1.5 — 2 Gyr. 

The high age of the secondary is interesting 
because it is inconsistent with age of the system 
as derived from stellar activity which is of order 
600 Myr (Melo et al. 2006). Recently, Schroter 
et al. (2011) reported on Chandra observations of 
an planetary transit of Corot 2A finding a simi- 
lar case. While they do not detect the transit in 
X-rays they find that the primary is X-ray bright 
with a luminosity ^ 1.9 • 10^^ erg s~^, indicating 
an age < 300 Myr. Meanwhile, a potential stel- 
lar companion was undetected down to a limit of 
Lx ~ 9 • 10^^ ergs s~^ which is inconsistent with 
the 300 Myr age and the distance of 270 pc. 

The beginning of the flare observed in 2009 in 
HD 189733 is at phase </> ^ 0.54, which coincides 
with a location 77° forward of the sub-planetary 
point and emerging to the day side of the star. 
This is almost exactly the location of the mag- 
netic sub-planetary point as calculated by Lanza 
(2008). The flare could be associated with the 
emergence of the foot-point of the magnetic col- 
umn to the earth facing side or a complex active 
region induced by magnetic SPI. Overall, the flare 
is associated with a change of the mean plasma 
temperature from ~ 0.5 keV to ^ 0.8 keV. The 
O VII triplet is in excess with respect to the best 
fit model. In the RGS spectra we observed that 
during the flare the inter-combination line seems 
to disappear, and the forbidden line is less lumi- 
nous. 

In this letter we report on the follow up obser- 
vation obtained in April 2011 at the same phase as 
in 2009, during an eclipse of the planet. Sections 
2 and 3 describe the observation and the results. 
Sect. 4 reports our conclusions. 

2. Observation and data analysis 

The observational setup mimics our previous 
observation made on May 18th 2009. We ob- 
tained a X-ray observation with XMM-Newton 
around the eclipse of the planet of HD 189733 
on April 30th 2011, starting at 23:14.20 (Obs Id: 
0672390201), and for a total duration of ~39.1 
ks. As in the previous observation, we used the 
Medium filter. The time between observations of 
HD 189733 was almost 2 years or exactly 61530.1 
ks (mid-eclipse to mid-eclipse). The period of the 
star is 11.953±0.009 days (Henry & Winn 2008), 



hence the star had rotated through 59.5795±0.045 
periods, while the planet had orbited it 321 times. 

For the reduction of the data we followed the 
same procedure as in Paper I, by using SAS ver. 
11.0 to extract events, light curves and spectra of 
HD 189733 recorded with EPIC camera and RGS. 
For fitting the spectra, we used 2-T VAPEC mod- 
els for the pre-flare phases (see description of light 
curve in Sect. 3.1) with different temperatures 
but coupled to have the same abundances. Abun- 
dances of Fe, Ne, and O were left free to vary 
while all other abundances were kept fixed at the 
solar values. For flare and post-flare phases we 
added a third VAPEC component keeping frozen 
the parameters of the first two VAPEC compo- 
nents. The abundances of the third component 
were linked and fixed. We obtained estimates of 
the temperature and emission measure of the flar- 
ing plasma, also following its fading after the flare. 
We applied the same procedure of best fit to the 
spectra obtained in 2007 (for the whole observa- 
tion), and 2009 (split in pre-flare, flare and post- 
flare phases) in order to compare the results. 

We merged the RGS spectra in the non-flare 
and in the flare state of all observations. Flare 
and non-flare state were fitted independently using 
Gaussian lines with a narrow intrinsic line width in 
the SHERPA fitting tool (Refsdal et al. 2011), so 
that the total line width is dominated by instru- 
mental broadening. Lines arc adjusted in wave- 
length to account for possible errors in the zero- 
point of the wavelength calibration, but the wave- 
length difference in multiplets is hold constant. 
Within the errors, all wavelength are compatible 
with the theoretical values. Instrumental back- 
ground and source continuum are assumed to be 
constant over a small wavelength region around 
the fitted lines. The fits are done using a Cash 
statistic, which takes into account the Poisson dis- 
tribution of counts, but the estimates of errors 
might be uncertain for very low count numbers. 

3. Results 

3.1. PN light curve cind spectra 

Fig. 1 shows the hght curve of EPIC PN in 
0.3-1.5 keV ^. The overall rate in 2011 is similar 



^Ephcmcridcs in http://var2.astr,o.cz/ETD/etd.php?STARNAME= 
HD189733&PLANET=B (Poddany et al. 2010) 
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Table 1: Values of best fit modeling to spectra before, during and after the flare, respectively for 2011 and 
2009 observations. Wc report also the fit of 2007 observation with two 2-T VAPEC components. The 
abundances of Fe, Ne and O in 2011 are derived from the pre flare phase, with values: Fe = 0.571q;^5, O = 
O.Sl^Q 07, Ne = 0.3^0^3^) respectively. The values of emission measures are given in units of 10^° cm~^. 
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to the PN rate recorded in 2009 and about twice 

the PN rate in 2007 observation (quiescent rate 
in 2009 and 2011 ~ 100 ± 12 ct ks"\ quiescent 
rate in 2007: ~ 60 ± 10 ct ks"^). The average 
energy of the spectrum in 2011 before and during 
the eclipse is 675 ± 20 eV, it has been 700 ± 10 eV 
in 2009 before the eclipse and 660 ± 10 eV during 
the eclipse. During the planetary eclipse we do 
not see the softening as in 2009. 

The most striking feature is the flare after the 
end of the planetary eclipse, which is analogous 
to the main flare seen in the 2009 observation. In 
2011, the flare starts at phase 0.52, while it starts 
at 0.54 in 2009. The duration of the flare (- 7 
ks, evaluated by eye), the peak rate (~ 230 ct 
ks~^ vs. ~ 210 ks~^ in 2009) and the presence 
of secondary impulses during the decay arc quite 
similar to those observed in 2009 as well. The 
detection of two flares within 120 ks (i.e. the sum 
of exposure times in 2007, 2009 and 2011 obser- 
vations) is interesting on its own, given that for 
active stars the typical rate of occurrence of bright 
flares is one every 500 ks (Wolk et al. 2005, and 
Caramazza et al. 2007) . 

Table 1 reports the best flt values of the spectra 
before, during and after this flare. For comparison, 
we report also the values obtained with the same 
scheme of best flt for 2007 and 2009 observations. 
Fig. 2 shows the PN spectra accumulated before, 
during, and after the bright flare observed in 2011. 

The pre-flare phase has a best flt with two ther- 
mal components, at 0.24 keV and 0.73 keV, re- 
spectively. In 2009 the values are lower (0.18 keV 



and 0.47 keV, respectively) and thus the corona 

appears colder than in 2011. 

During the flare the plasma has a marginally 
higher temperature, around 0.9 keV (flare tem- 
perature in 2009 was ^ 1 kcV). In facts, the flare 
spectrum shows an excess around 0.9 keV (Fig. 
2). 

In the post-flare phase the plasma cools essen- 
tially to the pre-flare temperature. The third com- 
ponent has a temperature similar to the hot com- 
ponent of the pre-flare and it is almost indistin- 
guishable from the latter. It is worth noting that 
the luminosity of the star remains slightly higher 
after the flare. This is observed in both 2011 and 
2009. 

As suggested in Paper I, this flare could arise 
from an active region which sits in a well deflned 
location on the stellar surface. It is plausible that 
in this region a strong magnetic field is present 
given its configuration. The periodic passage of 
the planet could trigger magnetic reconnections 
and thus strong fiares that could not arise else- 
where on the stellar surface. Fares et al. (2010) 
have published maps of the configuration of the 
magnetic field of HD 189733 in 2006, 2007 and 
2008. The field has a complex configuration with 
a toroidal component up to 40G in 2008 (the clos- 
est to the 2009 observation and the present one) 
and it changed configuration between 2006 and 
2008. The dipolar component is strongly non- 
axisymmetric in 2008. Cohen et al. (2011) have 
shown through MHD simulations based on these 
maps that magnetic SPI is possible whenever the 
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planetary and stellar alfvenic surfaces intersect 
each other. Given the non-symmetric shape of the 
stellar magnetic field, the interaction of planetary 
and stellar fields are not occurring continuously 
but at some phase determined by the planetary 
motion and the stellar rotation. 

The emission levels in X-rays in 2009 and 201 f 
are quite similar. On the other hand, these are 
^ 80% greater than in 2007. The observations 
suggests a cyclic activity which is in agreement 
with the changes in the magnetic field derived by 
Fares et al. (2010). We have no information on 
the topology of the magnetic field during last X- 
ray observation. The hypothesis that periodic oc- 
currence of strong variability is triggered by the 
planetary motion needs further observations to be 
assessed. 

3.2. RGS spectroscopy 

Merging the RGS spectra of two observations 
of XMM-Newton at the secondary eclipse and the 
one at the planetary transit (2007, PI: Wheatley), 
we have obtained the first high resolution spec- 
troscopy in X-rays of HD 189733. Table 2 shows 
line strengths of Fe XVII and the He-like Ne IX, 
O VII and the H-like O VIII. This latter is the only 
line which is significantly detected by both RGS 
detectors, the line fluxes are compatible within the 
statistical errors. Due to CCD failures, O VII is 
visible only in RGS 1, while Ne IX is visible only 
in RGS 2. 

3.2.1. Temperature 

In the flares, the Ne X emission increases by 
about a factor of three in agreement with the in- 
creased temperature seen in the light curve of me- 
dian energy and spectra of PN. The total emission 
of the He-like triplets is similar in flare and non- 
flare intervals. 

The ratio of lines from different ionization 
stages can be used as a temperature diagnos- 
tic (Mewe 1991). Figure 3 shows the ratio of 
O VIII/0 VII emission compared to the total lu- 
minosity in these lines. Blues squares show main- 
sequence (MS) stars from the sample of Ness et al. 
(2004), red triangles are classical T Tauri stars 
(CTTs), which are accreting pre-main sequence 
stars (Robrade & Schmitt 2007; Giinther 2011). 
CTTs show an excess of cool emission, which is 
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Fig. 1. — Top panel: light curve of PN in 0.3 - 1.5 
keV. The 2009 light curve (dotted line) is shown in 
phase with the 2011 light curve. Phase is marked 
on the top axis. The time bins are 400 sec. Bot- 
tom panel: median of energy (PI) as a function of 
time. The error bar of < PI > curve is about 
12 — 15 eV. Vertical lines mark 1st contact to 4th 
contact. Light curves of < PI > are smoothed by 
taking the median of the sample of 200 events and 
varying the sample by adding five new events and 
removing the five oldest ones (c.f. Paper I). 
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Table 2: Lines with significant detections. Error ranges are la or missing if unconstrained. 
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Fig. 2. — Spectra of PN before (black), during 
(red) and after the flare (blue). Differences in the 
shape are visible in the flare spectrum especially 
around 0.9 keV. 



likely powered by accretion shocks. HD 189733 
is among the stars in the sample -with the lo-west 
X-ray luminosity and O VIII/0 VII ratio. Thus, 
it is marginally cooler than MS stars of compara- 
ble luminosity, but its corona is over-dense like 
PMS accreting stars. All CTTs in Fig. 3 are sig- 
nificantly brighter than HD 189773, but this is 
a selection bias because only the brightest CTTs 
can be observed -with X-ray gratings. 

3.2.2. Density 

The ratio of the forbidden line / and the inter- 
combination line i in the He-like triplets is density 
sensitive (Mewe 1991; Porquet & Dubau 2000). 
In the lo-w-density, limit f/i is 3.5 for O VII 
(logrie < 9) and 3.2 for Ne IX {logn^ < 10) 
according to the CHIANTI database (Dere et al. 
1997, 2009). We performed Monte-Carlo simula- 
tions to account for the non-gaussian errors in the 
distribution (Giinther & Schmitt 2009, appendix 
A). The Ne IX lines are compatible -with the lo-w 
to moderate density limit in the non-flare state 
(95% confidence limit f/i > 2.2, i.e. logrie < 11), 
in the flare spectrum the error is too large to 
constrain the density. The non-detection of the 
O VII i line in the flare is consistent -with mod- 
erate densities, the 95% confidence lo-wer limit is 
f/i — 0.5 (logrie < 11.2). In the non-flare state, 
the 95% confidence upper and lo-wer boundaries 
on the f/i ratio are 1.8 (logng > 10.5) and 0.6 
(log Tie < 11.1). Generally, the coronal emission 
in MS stars is in the lo-w-density limit (Ness et al. 
2004), only in a fe-w bright flares higher densities 
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have been seen (e.g. on Proxima Can, Giidel et al. 
2002, 2004). 

The total luminosity in the He-like triplets 
changes little between the flare and the non-flare 
state. This could indicate that they originate in a 
dense region which is not aflFected by flare heating. 
However, the tentative change in the //i ratio can 
be interpreted as a change in density, where the 
hotter plasma seen in the flare state has a lower 
density than in the non-flare state. 

4. Conclusions 

We have analyzed a XMM-Newton observation 
of HD 189733 at the secondary transit of the 
planet. We observed a flare with characteristics 
very similar to that observed in the previous ob- 
servation at the same phase. The recurrence of 
such flares is explained by the following scenario: 
an active region is present at the same location 
on the stellar surface in both observations. The 
magnetic interaction with the planet is inducing a 
flaring activity in this region. The occurrence de- 
pends on the configuration of the stellar magnetic 
field and its strength. Further observations at the 
same phase can prove or reject this hypothesis. 

High resolution spectroscopy in quiescent and 
fiaring state shows that the corona of HD 189733 
is marginally cooler than MS stars with same oxy- 
gen line luminosity. We find a marginal change in 
f /i ratio of O VII lines in the flare state, implying 
a lower density during flares with respect to quies- 
cent state. In summary, the corona of HD 189733 
stands apart with respect to other coronae of Main 
Sequence stars without hot Jupiters. It is cold like 
the solar corona but it is dense and ten times more 
luminous and active. This could be related to SPI 
effects but more observations are needed to clarify 
the influence of its hot Jupiter on the corona. 

The XMM-Newton guest investigator pro- 
gram supported IP and HMG through grant 
NNX09AP46G. SJW was supported by NASA 
contract NAS8-03060 to the Chandra Science 
Center. VK acknowledges CXC NASA contract 
NAS8-39073, and OC acknowledges NASA LW- 
STRT grant NNG05GM44G. 
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Fig. 3. — Ratio of O VIII/0 VII emission com- 
pared to the total oxygen luminosity. The mea- 
surements for HD 189733 are shown with errors. 
The fiare measure has larger error bars because of 
the lower count number. See text for details (color 
in electronic version). 



This 2-column preprint was prepared with the AAS lAT^rjX 
macros v5.2. 
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